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Abstract. - Goldstone modes in the amorphous solid state, resulting from the spontaneous 
breaking of translational symmetry due to random localisation of particles, are discussed. Start- 
ing from a microscopic model with quenched disorder, the broken symmetry is identified to be 
that of relative translations of the replicas. Goldstone excitations, corresponding to pure shear 
deformations, are constructed from long wavelength distortions of the order parameter. The 
elastic free energy is computed, and it is shown that Goldstone fluctuations destroy localisation 
in two spatial dimensions, yielding a two-dimensional amorphous solid state characterised by 
power-law correlations. 



The elastic properties of gels and vulcanised matter have received considerable attention 
in the literature. Phenomenological arguments pp as well as microscopic models [2] have been 
used to compute, e.g., the shear modulus. The structural properties of gels have been discussed 
more recently 3 , building on the work of Edwards and collaborators 0] . It has been shown 
that in the gel phase a finite fraction of particles is localised around random positions with 
finite thermal excursions. What has not been analysed so far, and what constitutes the main 
focus of this Letter, is the connection between the spontaneously broken translational symme- 
try due to particle localisation and shear deformations, which are identified with Goldstone 
excitations as a consequence of the spontaneous breaking of translational symmetry. 

We first discuss symmetry-related equilibrium states in the general context of amorphous 
solids. We construct low-energy long-wavelength Goldstone fluctuations, which correspond to 
almost- uniform translations of the localised particles. Particular care is taken of the resid- 
ual symmetry of the amorphous solid state, which remains statistically homogeneous. We 
then apply these general ideas to a particular type of amorphous solid — gels and vulcanised 
matter — for which a replica field theory is available This allows us to identify the lost 
symmetry to be that of relative translations of the replicas, whilst common translations of 
all replicas remain as intact symmetries in the amorphous solid state. Goldstone fluctuations 
are described as long-wavelength distortions of the order parameter, which are reminiscent of 
the ripples of an interface between coexisting liquid and gas phases. The free energy of these 
fluctuations is computed within a Gaussian approximation, and identified with the elastic free 
energy of a macroscopically homogeneous isotropic amorphous solid. 
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Symmetry breaking and Goldstone modes. - We consider a system of iV distinguishable 
particles in a volume V interacting via potentials that are translationally invariant, e.g., 
additive pair potentials. In the amorphous solid phase a finite fraction of the particles are 
localised at random positions with thermal fluctuations of finite range. The system will be 
found in one of many possible equilibrium states. One such state is singled out arbitrarily, 



denoted by v and specified by nonzero local static density fluctuations (e 



k-R, 



Here, R,- is 



the position of particle j and (• • • )„ denotes a thermal average restricted to equilibrium state 
v. The amorphous state has macroscopic translational invariance, so that 



lim 



TV 
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?k-R, 



= 0, 
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where [• ■ • ] denotes a disorder average. Localisation in the equilibrium state v is detected by 
higher-order moments, the simplest one being the second: 



na(ki,k 2 ) = 



1 



^iki R, 



3 ik 2 -Rj\ 



= <5k!+k 2 ,o w(ki). 



(2) 



For this to be nonzero, macroscopic translational invariance requires ki +k 2 = 0, and w(ki) is 
a real-valued function. The second moment vanishes in the isotropic fluid phase and diagnoses 
the localisation of particles. Hence the second moment (and in fact all higher moments, as 
will be discussed below) serves as an order parameter for the transition from the fluid to the 
amorphous solid state. 

In a different equilibrium state, denoted by /i, all particles are translated with respect to 
the equilibrium state v by a displacement a' 1 , so that 



ik-R, 



ik-R, 



(3) 



Due to the translational invariance of the Hamiltonian there is in fact a manifold of symmetry- 
related states, generated by spatially uniform translations. How do we have to generalise 
the representation of the order parameter (J2J), in order to take into account the underlying 
translational symmetry of the Hamiltonian? For two equilibrium states, \i and p, the second 
moment is given by 



^' P (ki,k 2 ) = 



-y N 

N ^3 = 1 
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(4) 



We introduce a longitudinal wavevector ky = (ki 
(ki — k 2 )/2 and rewrite the above expression 



k 2 )/2 and a transverse wavevector kj 



ii£"'(ki,ka 



<Sk,,,o ^(kjj 



(5) 



in terms of a transverse displacement u^' p = a M — a p . Thereby, the second moment is repre- 
sented by its "magnitude" w(kj_) and the generator of transverse translations e lki ' u i P . 

The interpretation of the order parameter is more intuitive for the second moment as a 
function of real-space variables: V^2 ,p (xi, x 2 ) = cD(xi — x 2 — u P:P ). It does not depend on 
X|| = (xi +x 2 )/2, and hence can be visualised as a rectilinear hill, whose position is specified 
by a pair of pure states. In a simple model of particle localisation in D dimensions, a fraction 
Q of the particles are harmonically bound to random positions with random mean square 
displacements D £? , giving rise to 



w(kj 



q / dev(e)e- 



(6) 
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Fig. 1 - Second moment in real space: a hill with 
its ridge aligned in the xi + X2 direction. The 
thickness of the lines indicates the magnitude of 
the order parameter; the position of the ridge is 
determined by two equilibrium states. 
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Fig. 2 - Goldstone-distorted state in real space: 
the hill is displaced perpendicularly to the ridge 
by an amount u± that varies with position x\\ 
along the ridge. Note that the scale of this figure 
is much larger than that used for Fig. 



Here "P(£ 2 ) = [iV^ X^ ioc^(£ 2 — ^ s ^ e disorder-averaged distribution of squared locali- 
sation lengths of the iVi oc localised particles. The width of the hill is thus given by the typical 
value £ t yp of the localisation length. The structure in one space dimension is sketched in 

Fig. in 

This pattern of symmetry breaking suggests that the Goldstone excitations of an equi- 
librium state are constructed from it via xm- dependent translations of the hill in the xj_ = 
(xi — X2)/2 direction, i.e., ripples of the hill and its ridge (see Fig.[5J|. In Fourier space, the 
generalisation of Eq. (jjjj to inhomogeneous displacements reads 

Q(k 1 ,k 2 ) = y J dx|| e * r x ll+l k_ L .u i (x| ! ) (7) 
The corresponding representation in real space is completely equivalent: 

0(xi,X 2 ) = — ZD (x_L- U_l(X||)) . (8) 

Goldstone fluctuations in replica space. - In the following, we apply the general concepts 
of Goldstone excitations to a particular amorphous solid state — gels and vulcanised matter. 
For these systems a semi-microscopic model provides a framework for analysing the transition 
from the viscous fluid to the gel, as well as the structural and thermal properties of the gel, 
quantitatively. Here we use it to compute the free energy of Goldstone fluctuations, identify 
it with the elastic free energy, and thereby extract the shear modulus. 

The transition from the viscous fluid to the gel is controlled by the number of crosslinks, 
connecting randomly chosen pairs of monomers. These monomers can be part of a larger 
prefabricated structure, such as chains consisting of AI monomers or stars or dimers etc. For 
the following, the internal structure of the molecules is almost irrelevant; it will only determine 
certain parameter values. The average over the quenched disorder {i.e. the number and 
location of the crosslinks) is performed with help of the replica trick. The resulting, effectively 
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uniform, theory is formulated in terms of a field £l(x) defined over (l-fn)-fold replicated space 
x, so that the argument x means the collection of (1 + n) position D-vectors {x°, x 1 , . . . , x™}. 
According to the replica scheme of Deam and Edwards [3], the disorder- averaged free energy 
F is given by —FjT = lim n _>o (Zx+ n — Zi)/nZ\. where T is the temperature. The replica 
partition function Z\+ n is given as a functional integral Z\+ n ~ J T>Vl e~ Sn over the field 
fi(fc), where k = {k , k , . . . , k"} denotes wavevectors conjugate to x. 
The expectation value of the field is given by 

(**>>- (^E^iC,^),,- <9 > 

For this to be nonzero, macroscopic translational invariance requires X)q=o ^« = ®- The case 
with only one of the (n + 1) wavevectors nonzero corresponds to the (particle number) den- 
sity (Q. Density fluctuations are suppressed by excluded volume interactions, which tend to 
maintain homogeneity. Therefore the theory, in its simplest version, contains only fluctuating 
fields for which at least two of the D-component entries in the argument k = {k°, k , . . . , k" } 
are nonzero. We refer to such D(n + 1) dimensional vectors as the higher-order replica sector 
(HRS). Fields with I non- vanishing D-dimensional wavevectors correspond to the Ith moment 
of the local density, so in particular I = 2 yields the second moment discussed above. We 
consider a finite-dimensional system with a small number of nearest neighbours, implying 
that fluctuations of the local density are in general non-Gaussian and all moments need to be 
taken into account. 

The translational invariance of the system is reflected in the invariance of the effective 
Hamiltonian 6>n under independent translations of the replicas. In the amorphous solid state 
this symmetry is broken down to invariance under the subgroup of common translations of all 
the replicas. The manifold of symmetry-related classical solutions (i.e. saddle points of Sn) 
is given by 

fici(fc) =e ik ^S kh o u(k±), (10) 

in analogy to Eq. (J5J . We have decomposed wavevectors into longitudinal and transverse com- 
ponents with help of a complete orthonormal basis set in replica space {e Q }™ =0 . The "longi- 
tudinal" direction in (n + 1) dimensional replica space is the unit vector e = X)q=o ££> /vl+n, 
which is constant in replica space. All D(n + 1) dimensional vectors can be expanded in 
this basis: e.g., k — J^_ k Q e a and decomposed into longitudinal (||) and transverse (_L) 
components: 

k|| = _ Q k a /VTT^ and fc ± = fc-k||e. (11) 

Analogous expressions hold for x and u. Invariance under common translations of all replicas 
is guaranteed in the classical solution 11U|) by the Kronecker delta, which enforces the longi- 
tudinal wavevector k|| to be zero. The manifold of symmetry-related states is generated by 
purely transverse displacements. 

The order parameter (|1U|) exhibits the structure in replicated x-space discussed above 
for the second moment of the local density in different equilibrium states. The classical state 
f2 c i(x) can be visualised as a rectilinear hill in x-space with contours of constant height oriented 
parallel to xy e. Symmetry-related classical states are obtained by translating the hill rigidly, 
perpendicular to the ridge-line. Goldstone excitations from the classical state are constructed 
via X|| -dependent translations of the hill in the £_i_-direction [7], i.e., ripples of the hill and its 
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ridge: 

Vn(fe) = fdx {] e' ik n- x ii+^-^( x ii) w (A: ± ), (12a) 
Vtt{x) = D(xx-u_l(x||)). (12b) 

We are interested here in the low energy excitations only and, hence, consider Goldstone 
fluctuations, whose energy is expected to vanish with a power of the inverse wavelength of the 
ripples. We neglect all other fluctuations, which are expected to be massive. For consistency 
we thus should restrict ourselves to Goldstone fluctuations with wavelengths that are larger 
than a short distance cutoff l < . In particular we require that the displacement uj_(xii) varies 
on length scales larger than the width of the interface and take Z< = £typ- The physical content 
of this restriction is that, to be meaningful, elastic deformations should occur on length-scales 
longer than the typical localisation length of the particles in the elastic medium. 

As discussed after Eq. the critical fields of the theory reside in the HRS. Hence we 
want to make sure that the distortions of the order parameter remain in this sector. This is 
guaranteed if the distortions fulfill the incompressibility constraint: | det (<5 a ,& + d a u£) | = 1, 
which, for small amplitude distortions, reduces to d a u" — 0. 

Elastic free energy and shear modulus. - For the present purposes the essential term in 
the effective Hamiltonian 6>q is the "quadratic gradient" term, 

s Q = ^Vc $k 2 n(k)n(-k) + o (k A n 2 ,n 3 ,k 2 n 3 ) , (13) 

/cGHRS 

where c and £o are, respectively, the number-density and linear size of the underlying entities 
being constrained and HRS indicates that only wave vectors in the higher-order replica sector 
are to be included in the summation. Although Eq. (|13l) was originally obtained via an 
expansion valid close to the critical point, such a quadratic gradient term is expected to 
feature in a description valid deep inside the solid state. If one is seeking — as we are — the 
dominant contribution to the elastic free energy governing the displacement field u± [i.e. the 
leading-order term in u± and kn) then it is adequate to retain solely this quadratic gradient 
term. 

Inserting the fluctuating field into Sq and computing the increase, S u , due to the distortion 
111 (i.e. the elastic free energy) gives the following contribution, which arises solely from the 
quadratic "gradient" term in Sq: 

/n 
dx£ (a oU g'(x))(a.u?(x)) ) 

(1 + n) l+ ~ J nD 

Here, summation over repeated cartesian indices a, b is implied and the displacement vectors 
u Q have to fulfill the constraint ^ Q u Q (x) = 0. 

It was shown in Refs. |SJ|S] that in the classical state (|10fl 

/>oo 

w(*x) = 0/ deV(e)e~^ 2 , (15) 
Jo 

with the fraction of localised particles Q = 2ar/3g and the distribution of localisation lengths 
V(£, 2 ) — (£o/ ar £ 4 ) 7r (£o/ ar £ 2 ) i expressed in terms of a universal classical scaling function 



(14a) 
(14b) 
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7r(0). Here, ar and g are, respectively, the control parameter for the density of constraints 
and the nonlinear coupling constant. A semi- microscopic model of vulcanised macromolecular 
matter yields a — 1/2, r = (/i 2 — /j, 2 )//z 2 , £ 2 = M£p/2D and g — 1/6, where fi controls the 
mean number of constraints (and has critical value /x c = 1), £ P is the persistence length of the 
macromolecules, and M is the number of segments per macromolecule. If the network is built 
up from single monomers, corresponding to M — 1, then £ p is interpreted as the length of the 
tether connecting two point-like particles. 

We substitute this classical value of the order parameter into Eq. 1)1 4bl) and take the replica 
limit n — > to obtain the shear modulus 



The shear modulus increases continuously from zero with the third power of the distance from 
the gel point, in agreement with the previous microscopic calculation of Ref. 2 . 

A simple scaling argument is also consistent: The free energy density scales with the 
correlation length £ as T(~ D . This is to be compared with the elastic free energy density 
of Eq. 114a|) . Since d a Ub is dimensionless this implies po ~ C,~ D ~ t Du . At D — 6 we hnd 
agreement, which is to be expected since our calculation is based on a Gaussian expansion 
around the classical solution. 

Destruction of localisation in 2D and quasi long range order. - We discuss the effects 
of Goldstone fluctuations at the level of Gaussian fluctuations (|14a|) only. The elastic Green 
function is given by 



where k = k/|k|. Here, the integration over wave-numbers is restricted to 1/L < k/2-K < 1/Z< , 
where L is the linear size of the system. The mean value of the distorted classical state, 



(Vfi(fc))« / dx|| e lk »- x ii<e lfe -"^ x ii)) 5uW (fc ± ) = 5 k|| , o exp(-Tr D fci/2p ) u(k x ), (18) 



is determined by the local Green function Gab(0) = $a,b^D- In space dimension D > 3, 
is finite, so that the order parameter has a finite expectation value. The latter is reduced, 
due to Goldstone fluctuations, as one would expect: In addition to independent fluctuations 
of the localised particles there are also collective fluctuations — transverse phonons — which 
effectively increase the extent of spatial fluctuations of the localised particles. This is best 
seen by recalling the classical structure for uj, Eq. I|15fl : 



Order parameter correlations are given within the Gaussian approximation by 

(VSl(x t k x ) ra*(x',fcx)) w u 2 (k±) cxp (-?-(g ab ( ) - g ab ( x - x '))k ±a ■ k ±b ^j . (20) 

In D = 3, the asymptotics of the elastic Green function for large separation r = |x — x'| is 
well known, Q a b{^) ~ (^a,b + r a ff>)/87rr, implying for the correlator of the order parameter 




(16) 




(17) 




(19) 



(VSl(x,k±) VQ*{x', 



fcjj) w uj 2 (k±) exp 



( 



TT D k\ 



+ 



T{k\ + (kx ■ r) 2 ) 
8irp r 
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where r = r/r. 

In two space dimensions, Yr> increases logarithmically with the system size, so that the 
expectation value of the order parameter (|18fl vanishes in the macroscopic limit. Goldstone 
fluctuations destroy localisation in two space dimensions, in agreement with the Hohenberg- 
Mermin- Wagner theorem. Nevertheless, a two-dimensional amorphous solid persists, which is 
characterised by algebraic decay of order parameter correlations, i.e., quasi- long-range (ran- 
dom) order. The asymptotics of the clastic Green function for large spatial separation is given 
by Gab(0) — Gabi?) ~ 35^0,6 l n (r/Z<), implying an algebraic decay of the order-parameter 
correlations with a non-universal exponent r\^\ 

Tk 2 

(vQ( x ,fc x ) nr( x ',fc x )) ~ (r/uy^u 2 ^) with m = -±. (22) 

This scenario is similar to the theory of two-dimensional regular solids, taking into account 
harmonic phonons only, as discussed in Ref . [8] . The main difference is the absence of Bragg 
peaks at reciprocal lattice vectors. Instead, there is a divergence in the scattering function 
at k|| = 0, as can be seen by integrating the correlator, above, over all r = x — x' 
S(0,k x ) = f dr(VQ(r,k ± ) VQ*(0,k x )) ~ {L/l^ 2 -^. 

* * * 

Certain ideas concerning the nature of the vulcanization transition in two spatial dimen- 
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